H2 Line Emission 
Associated with the Formation of the First Stars 



Hiromi MIZUSAWA \ Ryoichi NISHI ^ and Kazuyuki OMUKAI ^ 

^Department of Physics, Niigata University, Ikarashi, Niigata 950-2181 
mizusawa@astro.sc.niigata-u.ac.jp 
nishi@astro.sc. niigata-u. ac.jp 
'^Division of Theoretical Astrophysics, National Astronomical Observatory, Mitaka, Tokyo 181-8588 

omukai@th.nao.ac.jp 

(Received ; accepted ) 
Abstract 

Molecular hydrogen line radiation emitted in formation events of first-generation 
stars are evaluated in a discussion of its detectability by future observational facilities. 
H2 luminosity evolution from the onset of prestellar collapse until the formation of a 
~ IOOM0 protostar is followed. Calculations are extended not only to the early phase 
of the runaway collapse but also to the later phase of accretion, whose observational 
features have not been studied before. Contrary to the runaway collapse phase, where 
the pure-rotational lines are always dominant, in the accretion phase rovibrational 
line emission becomes prominent. The maximum luminosity is also attained in the 
accretion phase for strong emission lines. The peak intensity of the strongest rovibra- 
tional line reaches ~ 10~^^W/m^, corresponding to the flux density of 10~^/iJy, for a 
source at the typical redshift of first-generation star formation, 1 -\- z = 20. Although 
the redshifted rovibrational H2 emission from such an epoch falls in the wavelength 
range of the next-generation infrared satellite. Space Infrared Telescope for Cosmology 
and Astrophysics, for exceeding the detection threshold 10^ such protostars are re- 
quired to reach the maximum luminosity simultaneously in a pregalactic cloud. It 
is improbable that this condition is satisfied in a realistic scenario of early structure 
formation. 

Key words: cosmology: early universe — galaxies: high-redshift — infrared: 
galaxies — stars: formation 
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1. Introduction 

First-generation stars can be a source of both the observed metals in the early universe 
and ultraviolet photons that cause the cosmological reionization. In fact, substantial metal 
enrichment has been observed in Lja clouds as early as z ~ 5 (Cowie et al. 1995; Songaila 
& Cowie 1996; Songaila 2001). Furthermore, the WMAP satellite has measured the electron 
scattering optical depth of Tg = 0.17 ±0.04, which corresponds to very early reionization at 
z ~ 17±5 (Bennett et al. 2003). The growing evidence, though still indirectly, suggest that 
first-generation stars did play a significant role in the early evolution of galaxies. 

Undisputedly the most decisive evidence for the first star formation is the direct detection 
of ongoing formation events. First star-forming clouds are considered to be bright in H2 line 
emission because the first stars are formed by gravitational collapse of primordial clouds, which 
is induced by H2 line cooling (Saslow & Zipoy 1967; Peebles & Dicke 1968; Matsuda, Sato, & 
Takeda 1969). 

In this paper, we evaluate the H2 line emission associated with first-generation star 
formation. Similar work includes that by Kamaya & Silk (2002) and by Ripamonti et al. 
(2002). They, however, have followed only the evolution in the earlier phase of star formation, 
namely the runaway collapse phase, and have not extended their calculations well into the later 
phase, namely the main accretion phase. Here, we study the observational appearance in both 
evolutionary phases. We have found that the maximum value of H2 line luminosity is indeed 
attained in the main accretion phase. 

Currently, several next-generation large observational facilities, including James Webb 
Space Telescope (JWST) ^, Space Infrared Telescope for Cosmology and Astrophysics (SPICA) ^, 
and Atacama Large Millimeter Array {ALMA) ^, are awaited. One of their major scientific goals 
is to explore the formation and evolution of pregalactic objects in the high-redshift universe. 
Among them, SPICA is the most relevant to our purpose since it has the best sensitivity in the 
mid- to far-infrared ranges, where the strongest H2 emission lines from the typical epoch for 
first star formation {z ~ 20) are redshifted. In this paper, we therefore discuss the detectability 
of those H2 line photons by SPICA. 

The outline of this paper is as follows: In section 2, the method of calculation is described. 
In section 3, the results for H2 emission from a forming first star are presented. Finally, in 
section 4, a brief summary is given and the detectability of first-generation stars by SPICA is 
discussed. 

^ http://www.jwst.nasa.gov/ 

^ http://www.ir.isas.ac.jp/SPICA/ 

^ http: / / www.alma.nrao.edu/ 
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2. The Model 



2.1. Density and Velocity Evolution 

From the viewpoint of dynamical evolution, it is convenient to consider the star formation 
process as two distinct phases, i.e., before and after the birth of the central protostar. These 
two phases are respectively called the runaway collapse phase and the mass accretion phase. 
Because of their vastly different dynamical natures, we treat these two phases separately and 
use different models for describing dynamics in each phase. 

2.1.1. The Collapse Phase 

Since the typical timescale for gravitational collapse, the free-fall time, depends inversely 
on density, the denser region collapses faster and becomes even denser relative to the rest. As a 
result, only the density in the central densest region increases, leaving the surrounding matter, 
or the envelope, almost unevolved. This central region has a fiat density distribution with a 



diameter given approximately by the central local Jeans length Aj = yc^n/Gp. On the other 
hand, the envelope has the power-law density distribution, p oc r"^/'-^"'^-' where r is the distance 
from the center and 7 = dlnp/dlnp is the effective adiabatic coefficient (Larson 1969). During 
the runaway collapse phase, the Jeans mass Mj = pXj^ decreases monotonically. Note, therefore, 
that matter that is a part of the central region but is off-set from the center at the beginning 
will be left behind at some instance during the collapse, and will be incorporated thereafter 
into the envelope. 

Considering this nature of dynamics, we use the following assumption for evolution in 
this phase: only when the matter is inside the central region, i.e., r < -Rcore, the physical 
quantities are assumed to evolve and the physical quantities outside do not evolve until the end 
of this phase. The radius of the central region is given by -Rcorc = where a is a correction 

factor of order unity and its value will be given below. In this way, only with the evolution of 
central quantities, we can Icnow the density and velocity structure of the entire cloud. 

The central evolution is assumed to be described by the free-fall relation modified by 
the pressure gradient: 

dp^_p_ . 
dt ^ ' 

where p is the density at the center and the free-fall time, tg, is 



Here, the modification due to the finite pressure gradient force is represented by the correction 
factor (3. 

Next we consider appropriate values for the correction factors a, and (3. Omukai & Nishi 
(1998; ON1998 hereafter) have studied the formation of the first stars by radiative hydrody- 
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namics assuming the spherical symmetry and have found that Larson-Penston-type similarity 
solution for 7 = 1.09 (Yahil 1983) is a good approximation for the actual collapse dynamics. 
For the Larson-Penston solutions, 

a = (3) 

where 5 is the asymptotic central value of the ratio of the pressure gradient force (per unit mass), 
Fp, to the gravitational force, Fg for the Larson-Penston-type similarity solution. The relations 
above can be derived from the following consideration: The ratio 6 is inversely proportional to 

i = fp/f,:.(-^)/(^)«^-^-«(-^r«a- (5) 



"'Core 



where P is the pressure and M^ore is the mass in the core. If -Rcore = Aj/2, i.e, the diameter of 
the core is equal to the Jeans length, the pressure force exactly balances the gravity, i.e., a = l. 
Thus 6 = If the pressure gradient force is neglected, the collapse timescale is the free-fall 



time tg = y37r/32Gp. In our case, the acceleration is reduced by a factor of 1 — 5 from the 
gravitational acceleration as a result of the pressure force. Then, by replacing the gravitational 
constant G in the free-fall timescale to (1 — 6)G, the collapse timescale is (1 — 6)^^^^tf[. That 
is, /3= (l-5)-^/2_ For 7 = 1.09, 

pressure force 
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0.78. (6) 



gravity force 
Then we use a = 1.13 and /? = 2.13. 

2.1.2. The Main Accretion Phase 

At very high density [n > 10^°cm~^), because of the lack of an efficient cooling mecha- 
nism, the temperature rises adiabatically and the increased pressure gradient force eventually 
halts the central collapse. At this stage, a hydrostatic core, or a protostar, is formed. The for- 
mation of the protostar marks the onset of the main accretion phase. Although the protostar is 
small (< 1O~^M0) at its birth, it grows by many orders of magnitude in mass due to accretion 
of the envelope matter. 

After protostar formation, the pressure gradient force in the envelope becomes negligible 
in comparison with the gravitational force. Then we assume that the infall of the envelope 
matter proceeds according to the free-fall rate after the protostar formation. The evolution of 
each mass shell is given by 
dv _ GMr 

Itt ~ ^ ' 

dr , , 



where r, f , and p are the distance from the center, the velocity, the density of a free-faUing 
mass shell, and is the mass contained inside the mass shell, respectively. In this case, the 
density of a mass shell increases as p oc r~^/^. We calculate the evolution of each mass shell until 
the beginning of H2 dissociation: because of the increasingly short timescale and reduced H2 
abundance, time-averaged H2 emissivity during the H2 dissociation is only a negligible fraction 
of the entire emissivity. 

2.2. Thermal and Chemical Evolution 

The thermal evolution is followed by applying the energy equation to each fluid element: 

(9) 




(10) 

7ad - 1 yUmH 

v= , (11) 

/iTTlH 

where e is the internal energy per unit mass, /i is the mean molecular weight and mn is the 
mass of the hydrogen nucleus. The net cooling rate Acooi is 

-^cool ~ -'^linc ~l~ ^cont ~(~ ■'^chcm; (1^) 

where Aune is the rate of cooling by H2 line emission, Acont is the net rate of cooling by continuum 
emission and Achem is the net rate of cooling by chemical reactions. 

The line cooling rate Aiinc is treated as in Omukai(2000), where the cooling rate owing 
to optically thick lines is reduced according to the escape probability formalism. The optical 
depth r of a line is evaluated locally by that for the shielding length 

4h = min(AS'th,«Aj/2); (13) 

At^D I 3At;D/5% (the runaway collapse phase) 

(dv/dr) 1 ^/2AvDr^^'^/y/GM (the main accretion phase), 

where Av^ is the velocity dispersion. As a source of the velocity dispersion, we consider only 



the thermal motion:, and then Avd = J2kT/ fimu. We use the relation for the homogeneous 



collapse V = r/{3f3ts) for the runaway collapse phase and the free-fall velocity v = -y/2GMr/r 
for the main accretion phase. 

The net continuum cooling rate Acont consists of the sum of cooling term by contin- 
uum emission and the heating term by continuum absorption. Both terms are treated as in 
Omukai(2001; eq. B27), but, for the heating term, we consider the radiation field from the 
central protostar. As mentioned below, because of the low effective temperature of the pho- 
tosphere, Acont contributes to cooling, rather than to heating, as a whole. In our calculation, 
we assume that the protostellar radiation is the diluted black body of 6000K from the pho- 
tospheric radius of lOOi?©. The photospheric radius is fixed at lOOi?© for simplicity. But, in 
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reality, it varies dramatically during the accretion phase, as pointed out by Omukai & Palla 
(2003). However, the radiation field of the photosphere infiuences the thermal state in falling 
mass shells only moderately because of the low effective temperature of the photosphere, and 
this simple treatment does not affect our results significantly. 

Major chemical reactions between hydrogen compounds H, H2, H+, H~, and 132" are 
solved. Included are reactions 1, 2, and 7-22 in Omukai (2000) with the same rate coefficients. 
Helium compounds are not included since they are thermally inactive in the temperature range 
in our calculation. The chemical cooling rate Achcm, which is mainly owing to H2 formation 
and dissociation, is also treated as in Omukai (2000). 

We neglect the external radiation field for simplicity. In the accretion phase, the radiation 
comes from the luminous central star. However, in the spherical accretion, the UV radiation 
is not strong, because of the low effective temperature of the star. Also the H2 column density 
well exceeds the value needed for the self-shielding of dissociating photons. Thus, in our model, 
photoionization and photodissociation are not important and not included. 

We use an assumption similar to that in section 2.1.1: the physical quantities, such as 
temperature and chemical composition, of the envelope remain constant until the end of the 
collapse phase. With this prescription, we can know the temperature and chemical composition 
distribution of the entire cloud as well as the density and velocity structure. Therefore, we can 
infer H2-line emission from a certain position in the envelope using the same method described 
above for calculating the emission from the central region. Summing the emissivity from the 
central region and from the envelope, we obtain the the H2 line luminosity of the entire star- 
forming core at each step. 

2.3. Initial Condition 

The fragmentation of metal-free gas clouds and formation of star-forming cores have 
been investigated by some authors (e.g., Uehara et al. 1996; Nakamura & Umemura 2002; 
Abel et al. 2002; Bromm et al. 2002). As the initial condition for our calculation, we adopt 
the typical values for the star-forming cores from Bromm et al.(2002): the number density 
riH = lO'^cm"^, temperature T = 200K, electron abundance y{e~~) = 10~^ and H2 abundance 
y{B.2) = 10~^. This condition corresponds to the Jeans mass of ~ lO^M©. 

3. H2 Luminosity from a Single Star-forming Core 

In Figure 1 the evolution of (a) the temperature and (b) the H2 concentration are shown 
against the number density. The evolution of the central region in the collapse phase is shown 
by the solid curve. For the accretion phase, the evolution of three representative mass shells 
is indicated. First, we describe evolution in the runaway collapse phase: In this phase, while 
the number density increases by about ten orders of magnitude, the temperature rises only by 
an order, owing to the radiative cooling by H2 (see Fig 1 a; e.g., Palla, Salpeter and Stabler 
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1983; ON1998). The rapid increase in the abundance of H2 around tt-h = lO^cm"^ is owing to 
the three-body reaction (Palla et al. 1983). As a resuh of this, most of the hydrogen in the 
central region has been in the molecular form before the number density reaches ~ lO^^cm"^. 
Next, let us consider the accretion phase: In this phase, the temperature is higher and the 
H2 concentration is lower than those in the runaway collapse phase for the same density. The 
reason for this is as follows: During the accretion phase, the compression timescale of each mass 
shell is shorter than that in the runaway collapse phase for the same density. This causes a 
rapid temperature rise. As a result, the H2 dissociation starts before the formation of sufficient 
molecular hydrogen for cooling. 

The four strongest H2 lines in our calculation are rovibrational lines of l-OQ(l), 1-00(3), 
pure rotational lines of 0-0S(4) and 0-0S(3)^, whose wavelengths at rest are 2.34/im, 2.69/im, 
8.27/im, 10.03/im, respectively. The evolution of those line luminosities is shown in Figure 
2. During the runaway collapse phase, pure-rotational lines are stronger than rovibrational 
lines. In contrast to this, in the accretion phase, the rovibrational lines become stronger than 
pure-rotational lines because of the high temperature and density in this phase. Note that the 
luminosities of both rovibrational and pure-rotational lines reach the maximum in the accretion 
phase. The peak values of the line luminosities Lpeak and the mass of the protostar at this epoch 
Mpeak are shown in Table I. 



Table I. maximum line luminosities 



line 


A(/im) 


W(10''erg/s) 


Mpeak (Mq) 


l-OQ(l) 


2.34 


1.73 


12.9 


1-00(3) 


2.69 


1.71 


11.2 


0-0S(4) 


8.27 


0.78 


7.4 


0-0S(3) 


10.03 


1.05 


10.5 



Here we discuss why the line luminosities are stronger in the accretion phase than in 
the collapse phase. The possible sources of energy radiated by the H2 emission include the 
gravitational energy liberated by the infall and the chemical binding energy liberated by H2 
formation. Before the peak of luminosity, the main source is the latter. The majority of the 
former turns into kinetic energy without being converted to radiation energy. On the other 
hand, most of the chemical binding energy liberated by H2 formation is converted into thermal 
energy on the spot, since the number density is high enough to assure the local thermodynamic 
equilibrium (n ^ lO'^cm"^) in our calculation. The evolution of the liberation rate of chemical 
energy and the H2 radiation emissivity for the mass shell that contains about I3M0 inside are 
shown in Figure 3. Hereafter we write this mass shell as "mass shell (I3M0)". As is seen in the 

^ For example, the meaning of l-OQ(l) is as follows. " 1-0 " denotes vibrational transition from u = 1 to f = 0, 
Q denotes that the change of rotational level, A J, is 0. Similarly, S and O denote AJ — —2, +2, respectively. 
The number in the bracket denotes rotational level after the transition. 
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figure, the majority of tlie radiated energy is indeed coming from the hberated chemical binding 
energy. We compare here the hberation rate of chemical energy in the collapse phase with that 
in the accretion phase. In Figure 4, we show the chemical energy liberated per unit time at four 
selected instances, the final stage of the collapse phase and three steps of the accretion phase, 
against the radial distance. The liberation rate of chemical energy is higher in the accretion 
phase than in the collapse phase. The reason for this is as follows: The radius of the H2 forming 
region, which is defined as the radius where the H2 formation rate is maximum, decreases from 
~ 3 X lO^^(cm) in the collapse phase to ~ 6 x lO^^(cm) when the mass of the central protostar is 
~ 5Mq in the accretion phase. At the same time the inflow velocity and thus the mass entering 
into the H2 formation region per unit time becomes larger in the accretion phase, i.e. the total 
number of H2 forming per unit time increases. Since the liberated chemical binding energy is 
converted into radiation efficiently, the luminosity in the accretion phase exceeds that in the 
collapse phase. 

The luminosities of rovibrational lines reach their maximum at a protostellar mass of 
about 13Mq (Figure 2). Here, we discuss why this luminosity peak appears. In Figure 5, we 
show the evolution of the mass shell (13Mq) inside, against the radial distance. For comparison, 
we also plot some of the other mass shells. When the mass of the protostar is small, despite the 
large emissivity per unit mass, the mass of the region appropriate for vibrational lines, where the 
temperature is between lOOOK and 2000K, is small. On the other hand, at a later phase, owing 
to the fast contraction due to the strong gravity of the central protostar, the temperature rapidly 
reaches the H2 dissociation value 2000K) before sufficient H2 formation. The radiative 
emissivity is smaller, although the mass of the emission region is larger. Consequently, the 
luminosity peak appears at some epoch in between. 

In the outer envelope where density is low (~ 10^ — lO^^cm"^), the temperature of in- 
falling mass shells in the collapse phase is higher by a factor of ^ 2 than the value of the 
Larson-Penston-type similarity solution for 7 = 1.09. Therefore, the mass accretion rate, which 
depends on the temperature as tI, dM/dt cl/G (x T^/^ (e.g., Stahler et al. 1986), becomes 
larger than that derived from the similarity solution by a factor of ^ 4 after the protostellar 
mass exceeds ~ 3M0. The mass accretion rate so calculated decreases slowly from the initial 
value of ~ 0.05MQ/yr to ~ 0.02MQ/yr at the protostar mass of IOOMq. 

4. Discussion and Conclusion 

We have studied the evolution of H2 line emission from a forming first star, starting from 
a prestellar dense core up to the formation of a IOOMq protostar. After the formation of the 
central protostar, the H2 line luminosity increases and luminosities of strong lines reach their 
maximum when the central protostellar mass is ~ 13M0. The four strongest H2 lines are the 
rovibrational lines of l-OQ(l) (2.34/im), 1-00(3) (2.69/im), the pure-rotational lines, 0-0S(4) 
(8.27/im), and 0-0S(3) (10.03/im). During the accretion phase, vibrational lines of l-OQ(l) and 
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1-00(3) are the strongest, while in the collapse phase the strongest lines are the pure-rotational 
lines. 

Strong H2 line emission is expected also from diffuse regions in low-metallicity {Z < 
O.OIZ0) forming galaxies (e.g., Shchekinov & EnteF 1985; Shchekinov 1991; Omukai & 
Kitayama 2003). However, the dominant lines from those sources are pure-rotational lines be- 
cause of the lower density {n ~ 10'^~'^cm~^). On the other hand, in the late stage of the first star 
formation, H2 rovibrational lines are dominant, because of both high density (n ~ 10^°~^^cm~^) 
and temperature (T > lOOOK) in the H2 emitting region. Therefore, detection of H2 rovibra- 
tional lines can be a useful method to identify first-generation star forming clouds. 

Several large observational projects using next-generation facilities are now in progress, 
including JWST, SPICA and ALMA. The wavelength ranges of these three facilities are pre- 
sented in Table II. 



Table II. The wavelength ranges of future telescopes 



JWST 


SPICA 


ALMA 


0.6-28/im 


5 - 200/im 


300/im — 3mm 



For typical formation epochs of the first stars z ~ 20 (e.g., Abel et al. 2002), the emitted H2 
rovibrational line photons are redshifted to the far-infrared wavelength. As is seen in Table II, 
SPICA has the deepest detection limit in this range and thus is the most suitable to our purpose. 
SPICA is a next-generation infrared satellite, which is planned to have a cooled (4.5K) and 
large (3.5m) telescope and to be highly sensitive at the mid- and far-infrared wavelengths. The 
line detection limit of SPICA is about 10~^'^W/m^ in the far-infrared wavelength A~40/im with 
a spectroscopic resolution R = 3000 and observation time of one hour. Since a signal to noise 
ratio is proportional to ^/t, the sensitivity improves, of course, with longer integration time. 
Nonetheless, in the mid- and far-infrared wavelengths, the high background of the zodiacal light 
limits the sensitivity of SPICA to 10~^^^~^^-*W/m^ (private communication H. Matsuhara). 

The observed wavelengths of the four strongest lines redshifted from the emission at 
1 + z = 20, the typical epoch of first star formation, are given in Table III. 



Table III. The Observed Wavelengths of Strong H2 Lines 



line 


l-OQ(l) 


1-00(3) 


0-0S(4) 


0-0S(3) 


rest frame 


2.34/im 


2.69/im 


8.27/im 


10.03/im 


redshifted (1 + 2 = 20) 


46.8/im 


53.8/im 


165.4/im 


200.6/im 



The wavelengths of the redshifted pure-rotational lines shown above are outside the range of 
any future telescope given in Table III. On the other hand, the wavelengths of the rovibrational 
lines, which are intrinsically stronger than pure-rotational lines, fall in the range where the 
sensitivity of SPICA is excellent. 

By using the value of peak luminosity Lp^ak in Table I, the peak intensity of the strongest 
rovibrational lines from 1 + ^ = 20 is 
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peak 



1.73 X lO^'^erg/s 



0.32 X 10"2^W/m2 



(15) 



4vrZ^i+.=20 



47r(2.13 X 10iipc)2 



where -Di+2=2o = 2.13 x lO^^pc is the luminosity distance to 1 + 2 = 20 given by the formula 



where Hq is the Hubble parameter, with the values Hq = 70kms~^Mpc~ , = 0.7, f^M = 0.3. 
This peak intensity corresponds to the peak flux density of ~ 10~^/iJy using the line width 
vqv/c ~ lO^^Hz, where vq is the line-center frequency for the observer and v is the velocity of 
the mass shell. Note that in considering a cluster of star-forming cores, a larger value of the 
line width should be taken because relative velocity between the cores is usually larger than 
the thermal velocity in a core. In this case, the peak flux density per unit core will be lower 
than the value given above. 

Taking into account the limiting sensitivity 10~'-^^~^^''W/m^ and the maximum rovibra- 
tional line flux ~ 10"^^W/m^, we conclude that SPICA is able to observe a cluster of forming 
first-generation stars if it is contains more than 10^ sources with their luminosity near the peak 
value. 

Next, we consider whether or not this situation is possible in a realistic cosmological 
context. The epoch of the peak luminosity lasts only for Atpoak ~ lO'^yrs (see Fig 2), which is 
by far shorter than the typical free-fall time in the pregalactic clouds (say, ^ lO^yrs). Even in 
the extreme case of starburst, star formation is a prolonged event that continues at least a free- 
fall timescale of the pregalactic cloud. Thus, at a given time, at most only Atpcak/^s ~ 10~^ of 
the total gas mass of the pregalactic cloud is a total mass of star forming cores approximately at 
peak luminosity. To reach the luminosity peak, the first-star forming cores must be as massive 
as 13Mq (§3). The mass scale of actual cores is probably more massive, typically IO^-^Mq, 
according to numerical simulations of the fragmentation of primordial clouds (e.g., Bromm et 
al. 2002). Thus the total mass of 10'' cores is ~ lO^M©. Since only 10~^ of all the cores are 
at peak luminosity, in order for 10^ cores to reach peak luminosity at the same time, the total 
cloud mass must be as high as IO^^Mq. Note also that in this case the star formation rate is 
required to be an improbably high value of lO^MQ/yr. It is true that for pure-rotational lines, 
i.e., 0-0S(3), the epoch in which the luminosity is higher than lO^^erg/s continues for as long 
as lO^yrs. Although this long emission period reduces the required star formation rate, it is 
still improbably high (lO^MQ/yr). On the other hand, in order to emit radiation mainly in II2 
lines during the prestellar collapse, especially for the very high density region, the metallicity 
in the pregalactic clouds must be lower than ~ 10~^Zq (Omukai 2000). However, formation of 
Galactic scale objects with such low metallicity is probably excluded in the context of standard 
cosmology (Scannapieco, Schneider, & Ferrara 2003). In conclusion, it is improbable that 
detection of H2 line emission from forming first stars can be accomplished with SPICA. 

In this paper, the accretion flow has been assumed to be spherical. In reality, however, 




(16) 
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an accretion disk should form owing to an angular momentum barrier. In this case, the total 
energy emitted in H2 lines might increase, since the released gravitational energy is converted 
more efficiently into radiation in the accretion disk. At the same time, H2 is more susceptible 
to photodissociation owing to direct irradiation from the star. To estimate the H2 luminosity 
more precisely, it is necessary to study effects of the disk in the accretion phase. We leave 
this problem for future work. Recently, Bromm & Loeb(2004) performed a three-dimensional 
numerical simulation and derived the mass accretion rate onto the protostar. Their accretion 
rate is similar to our value until ~ 10^ yr after the protostar formation. After this time, 
their accretion rate drops more rapidly than ours. This difference is probably caused by more 
realistic evolution that they calculated, e.g. quasi-static contraction in the early evolution 
and/or the complex accretion flow owing to 3D effects etc. In our calculation the maximum 
values of the H2 luminosities are reached at ~ 300 yr after the protostar formation. At this 
epoch, our accretion rate does not differ much from their value. Therefore, the peak value of 
H2 luminosity is expected to be almost the same even if their accretion rate is adopted. After 
~ 10^ yr of protostar formation, we are probably overestimating the accretion rate and thus 
the H2 luminosity. Nevertheless, the H2 luminosity at this time is already weak and does not 
important for our discussion on observability by future facilities. 

Contrary to our claim, Kamaya & Silk (2002) and Ripamonti et al.(2002) have con- 
cluded that future observational facihties, e.g., JWST, ALMA, ASTRO-F, are able to detect 
H2 emission from the ffist star-formation sites. Here we stress that our estimated value of the 
peak luminosities of H2 lines is larger than their value, which is attained for the earlier phase. 
Moreover, in the mid- to far-infrared wavelengths, the sensitivity of SPICA that we use here 
to estimate the detectability exceeds that of the other facilities that these authors considered. 
The reason for the apparent contradiction is as follows: in converting the emitted H2 luminosity 
to the observed flux, they used the luminosity distance underestimated by a factor of 1 + z, 
where z is the redshift of the source. This error propagates to the final value of the observed 
flux, which is overestimated by (1 + zY . Note that, anyway, our value of the intrinsic luminos- 
ity in the runaway collapse phase is consistent with Kamaya & Silk (2002) and Ripamonti et 
al.(2002). 

Finally, we touch on another future telescope, the Single Aperture Far-Infrared 
Observatory (SAFIR) ^. Thanks to its larger mirror size (8 m), the line detection limit of 
SAFIR is better than that of SPICA and the background limit of ~ lO-^^W/m^ at A ~ 40/im 
can be reached in one hour with spectroscopy of i? = 1000. Recall that the detection limit 
we used in this paper is also the background limit ~ 10~^^W/m^, which is set by the zodiacal 
light. Therefore, even with SAFIR we conclude that detecting H2 line emission from forming 
ffist stars is difficult. 

The authors would like to thank Tetsu Kitayama, Hajime Susa and Masayuki Umemura 
^ http://safir.jpl.nasa.gov/whatls/index.asp 
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Fig. 1. (a) The temperature evolution and (b) the evolution of the H2 (bold Une) and H (thin hnc) 
concentration against the number density. For the runaway collapse phase, the evolutionary trajectory of 
the central part is drawn with a solid Une. For the accretion phase, the evolutionary trajectories of the 
mass shells of A//M0 = 1 (long-dashed lines), 13 (short-dashed lines) and 100(dotted lines) are shown. 
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Fig. 2. The luminosity evolution of the four strongest H2 lines, the rovibrational lines of l-OQ(l) (solid 
lines), and 1-00(3) (long-dashed lines), and the pure rotational lines of 0-0S(4) (short-dashed lines), and 
0-0S(3) (dotted lines). The horizontal axis is the time since the onset of the collapse. The evolution up 
to the formation of a lO^M© protostar is shown. The formation of the central protostar corresponds to 
t = 0. Note that all lines attain their maximum luminosity at t > 0, i.e., in the main accretion phase. Also 
the luminosities of rovibrational lines become stronger than that of pure-rotational ones in this phase. 
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Fig. 3. The evolution of the liberation rate of chemical binding energy (dotted line) and H2 line emissivity 
(solid line) for the mass shell (13Mq). The time t = is set when the mass shell (ISM©) is supposed to 
reach the center (r = 0). Because of the finite radius of the star, the mass shell reaches the stellar surface 
before t = Q. 
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Fig. 4. Liberation rate of chemical binding energy, at the final stage of the collapse phase (solid line), 
when the mass shells of Mi./Mq = 5 (long-dashed line), 13 (short-dashed hne) and 50 (dotted hue) fall on 
the center, respectively. The horizontal axis is the distance from the center. The area surrounded by each 
line and the horizontal axis indicates the chemical energy liberated per unit time. 
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Fig. 5. The evoluton of (a) number density, (b) temperature, (c) H2 concentration of selected mass shells 
against the radial distance from the central protostar in the main accretion phase. The values for the mass 
shell (I3M0) are shown by the thick solid line. Recall that the radii of the mass shells decrease with time 
since the envelope matter is falling onto the protostar. For comparison, we also plot the values of other 
two mass shells that contain IM0 (dashed line) and IOOMq (dotted line) inside. 
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